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Abstract 

In this work we present a theoretical and numerical study of the behaviour 
of the maximum Lyapunov exponent for a generic coupled-map-lattice in the 
weak-coupling regime. We explain the observed results by introducing a suit- 
able continuous-time formulation of the tangent dynamics. The first general 
result is that the deviation of the Lyapunov exponent from the uncoupled-limit 
limit is function of a single scaling parameter which, in the case of strictly 
positive multipliers, is the ratio of the coupling strength with the variance 
of local multipliers. Moreover, we find an approximate analytic expression 
for the Lyapunov exponent by mapping the problem onto the evolution of 
a chain of nonlinear Langevin equations, which are eventually reduced to a 
single stochastic equation. The probability distribution of this dynamical 
equation provides an excellent description for the behaviour of the Lyapunov 
exponent. Furthermore, multipliers with random signs are considered as well, 
finding that the Lyapunov exponent still depends on a single scaling param- 
eter, which however has a different expression. 

Short title: Maximum Lyapunov exponent in coupled-map lattices 
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I. INTRODUCTION 



Coupled-map lattices (CML) represent an interesting class of models for the investigation 
of several spatio-temporal phenomena, ranging from pattern formation to synchornization 
and to spatio-temporal chaos. Even though the discreteness of both space and time variables 
makes CML more amenable to numerical simulations than partial differential equations, the 
development of analytical techniques remains a difficult task. As usual, when dealing with 
hard problems, it is convenient to start from the identification of some relatively simple limit 
and thereby developing a perturbative approach. In the case of lattice dynamics, there are 
two opposite limits that are worth being investigated: the weak- and strong-coupling regime. 
In the former case, one can use all the knowledge acquired about low-dymensional systems 
to predict the dynamical properties when spatial directions are added. In this spirit, general 
theorems have been formulated and proved both about the structure of the invariant measure 
^ and about the existence of travelling localized excitations 0. The weak-coupling limit is 
interesting also in connections with the synchronization of chaotic attractors, a problem that 
can be effectively studied by looking at the behaviour of the Lyapunov exponents 0. In the 
opposite limit, one expects a slow spatial dependence and, correspondingly, the existence of 
a few, dynamically active, degrees of freedom. This is the spirit that has led to studying 
different truncations of partial differential equations. 

In this paper, we shall consider the most common indicator of chaos, the maximum 
Lyapunov exponent (MLE) in a lattice of coupled maps. It is known that a (weak) spatial 
coupling has a threefold consequence on the MLE: i) it naturally modifies the invariant 
measure; ii) it induces correlations among the local multipliers by coupling neighbouring 
sites; Hi) it modifies the dynamics of perturbations by inducing a coupling in tangent space 
as well. The last phenomenon is the most interesting one, since it leads to counterintuitive 
effects like an increase of the MLE even in the presence of a stabilizing coupling. 

In some previous papers 0,0, this problem has been addressed with a specific interest 
for the scaling behaviour of the MLE. However, the various approaches implemented so far 
have not been able to go beyond a qualitative explanation of the dependence on the coupling 
strength. Here, instead, we aim at presenting a fully quantitative, though approximate, 
treatment for the MLE in the small coupling regime. We restrict our analysis to the usual 
diffusive coupling scheme, but we are confident that the present approach can be effectively 
adapted to different (short-range) interaction schemes. More specifically, we refer to the 
dynamical equation 

Xi{t + 1) = f{y,{t + l)) 

y,{t + 1) = £x,_i(t) + (1 - 2e)xi{t) + exi+i{t). (1) 
The corresponding evolution equation in the tangent space is 

Ui{t + 1) = mi{t)^yeu^^l{t) + (1 - 2e)ui{t) + eui+i{t)^, (2) 

where mi{t) = f'{yi{t + 1)). From Eq. (||), we see that the last of the above mentioned 
three effects of the spatial coupling can be isolated and studied separately from the other 
ones. It is sufficient to assume that the distribution of multipliers is independent of the 
the coupling strength e and that the mi{t)^s are 5-correlated both in space and time. These 
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approximations are equivalent to introducing a random-matrix approximation in full analogy 
with Refs. 

The first author who investigated the effect of a small coupling on a chaotic dynamics was 
Daido , who studied numerically two coupled maps as well as two continuous-time chaotic 
oscillators. He also attempted a combined analytical and numerical study to explain the 
observed behaviour, without however being able to go beyond the prediction of the scaling 
behaviour. 

Later on, the problem of (infinitely many) coupled maps was considered in Ref. [Q, 
by exploiting the analogy with the statistical mechanics of directed polymers in random 
media. Indeed, Eq. (0) can be also read as the recursive equation for the partition function 
of a polymer of length t which grows by adding each new monomer no farther than 
one site from the last one. A possibly relevant difference between the two problems comes 
from the "Boltzmann weight" rrii which is necessarily positive in the polymer case (being 
a probability), while it can be negative in a CML (mj being the derivative of the map /). 
This is a first issue that makes the problem (^ more difficult to be treated and it is the 
reason why previous studies have been restricted to the case of strictly positive multipliers 
0J^. In fact, without entering the mathematical treatment, one can see that if rrii can be 
either positive or negative, Ui{t) is no longer positive definite and partial cancellations can 
occur in the iteration of the recursive relation. 

The efforts made in 0] to estimate the MLE consisted in developing a mean field analysis 
on the basis of the equivalence between the MLE in CML and the free-energy in directed 
polymers. By thus using the approaches developed in , it was possible to show that the 
spatial coupling induces an increase of the MLE from the "quenched" average Aq = (logmj) 
(corresponding to the absence of a spatial coupling) towards the "annealed" average Aq = 
In(mj), holding above some critical coupling value. While some features of this scenario were 
qualitatively confirmed by the numerical simulations (as, e.g., the increase of the MLE), no 
evidence of the phase transition was actually found. A perturbative technique, developed 
later on to improve the previous estimates has shown that the transition point slowly shifts 
towards larger coupling values, possibly disappearing in the asymptotic limit. Nevertheless, 
the extremely slow convergence of the estimates of the MLE to the values numerically 
observed, makes a general implementation of this approach not very appealing. Moreover, 
we should also recall that the analogy with directed polymers does not even allow an exact 
prediction of the scaling behaviour of the MLE, in so far as it indicates the existence of an 
additional, extremely weak, dependence on the coupling strength which seems to be absent 
in the outcome of direct numerical simulations. 

It is also worth recalling the analogy between the behaviour of the MLE and the evolution 
of rough interfaces. By interpreting the logarithm of (the amplitude of) the perturbation 
as the height of an abstract interface, the Lyapunov exponent becomes equivalent to the 
velocity of one such interface P,p]0| . However, this analogy is of no utility in the present case, 
since the deviation of the MLE from the uncoupled limit cannot be determined by studying 
the corresponding Kardar-Parisi-Zhang (KPZ) equation as already remarked in [Q. In a 
sense, the MLE corrections are connected to non-trivial deviations from a KPZ behaviour 
over short temporal and spatial scales. 

In the present paper we derive approximate but analytical expression for the MLE, by 
mapping the original problem onto that of a chain of continuous-time, nonlinear Langevin 
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equations. Such a set of equation is then reduced to a single stochastic equation whose 
solution yields an expression for the MLE that is in good agreement with direct numericals 
simulations. As the whole approach does not make use of the specific structure of the initial 
equations, we are rather confident that it can be repeated for other types of couplings, 
the only difference being presumably the structure of the deterministic force in the final 
stochastic equation. Moreover, we would like to point out that the mapping onto continuous- 
time equations indicates that the initial discreteness of the time variable is not a distinctive 
feature and we can imagine that a similar approach works also in the case of coupled chaotic 
oscillators. Our hope is reinforced by the observation that some equations obtained in the 
present framework can be derived also in the case of two weakly coupled differential equations 
0. 

The paper is organized as follows. In Sec. II, we briefly introduce the problem and 
some notations. In Sec. Ill we discuss the properties of tangent dynamics in the case of 
strictly positive multipliers: the analytical treatment is followed by a comparison with the 
numerical results. In Sect. IV we extend the analytical treatment to the case of multipliers 
with random signs. Finally, in Sect. V we summarize the main results and comment about 
the still open problems. The two appendices are devoted to the small-noise limit in the case 
of strictly positive multipliers and, respectively, to the small coupling regime in the general 
case of random signs. 



II. PRELIMINARY TREATMENT 



In this section, we formulate the problem of computing the MLE under the assumption 
of a small coupling strength. The first step consists in introducing the ratio between the 
perturbation amplitude in two consecutive site. 



which allows writing Eq. (0) as 

Ui{t+1) 



In 



Ui{t) 



In + In 



R,{t) 



(3) 



(4) 



where the average of the l.h.s. is nothing but the MLE A{e), while the r.h.s. is naturally 
expressed as the sume of the 0-coupling value plus the correction term induced by spatial 
interactions. 

A more convenient way of writing the Lyapunov exponent is obtained by transforming 
the smallness parameter as 



which leads to 



A{e) 



In 



Ui{t + 1 



Ui(t) 



A(0)+ln(l-2£) + (5A(7), 



(5) 



(6) 



where the Lyapunov correction is split into two parts, a multiplier-independent term and a 
non-trivial contribution 
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M(7) 



In 



1 + 7/?,+i + ^ 



(7) 



where (■) represents the time-average along the trajectory generated by Eq.(|l]) or, equiva- 
lently, the ensemble-average if ergodicity holds. 

Eq. (|^ tells us that the determination of the MLE requires the prior knowledge of the 
invariant measure of the stochastic process Ri{t) on two consecutive sites. By expanding 
the logarithm in powers of 7, the correction to the MLE can be expressed in terms of all the 
momenta of the probability distribution, 



(8) 



It is, therefore, necessary to formulate the dynamical equation in terms of the variable Ri. 
This can be easily done from Eqs. (|^,|) 



Riit + 1) = ixi{t) 



Ri{t)+-fRiit) Ri+^it)+j 



(9) 



l + 7i?i(t)+7/i?i„i(t) 
where we have introduced the stochastic process 

IJi{t)=mi{t)/mi^i{t) (10) 
whose geometric average is equal to 1, /i, being the ratio of two i.i.d. processes. 



III. POSITIVE MULTIPLIERS 

We first consider strictly positive multipliers, while the more general case is addressed in 
the following section. In fact, while the two cases require a somehow different treatment, a 
discussion of the former problem allows introducing several tools that turn out to be useful 
also in the latter context. 



A. Theory 

If the yUj's are positive definite, the ratios i?j's remain positive whenever initialized as 
such. This allows introducing the variable Wi = In Ri without the need of dealing with the 
sign of Ri- The introduction of Wi is convenient in that it transforms the original problem 
into a stochastic process with additive rather than multiplicative noise. With no restriction 
other than the positiveness of the multipliers, we obtain 

Wi{t + 1)~ Wiit) = \nfii{t) + In |l + 7e"'-+iW + ^e-"''^! - 

ln|l + 7e"''W +7e-"''-^(*)|. (11) 

In the limit of small 7, the arguments of the logarithms differ slightly from 1, provided that 
Wi does not deviate too much from 0. Under this approximation (that can be checked a pos- 
teriori), we can expand the logarithms, retaining the leading terms and thereby introducing 



5 



a continuous-time representation (because of the smallness of the deterministic variations). 
One obtains 

Wi = -27 sinh(wi) + 7 (^e"''+^ - e""''-^^ + ^i{t) (12) 

where = In/ij, has zero average and correlation function, 

imUt')) = 4a^'5(t - t')i6,, - ^5.±i,,) (13) 

where cx^ is the variance of the logarithm of the multiplier rrii and there is a factor 4 (instead 
of the usual 2) since the noise term is the "sum" of two i.i.d. processes {^i = Inmj — lnmj_i). 
Spatial correlations of the stochastic forces are all zero except those of neighbouring sites, 
a feature induced by the very definition of C,i as the sum of two processes in neighbouring 
sites. 

It is important to notice that the only parameter of the sequence of multipliers which 
eventually contributes to the correction of the MLE is precisely the amplitude of the fluc- 
tuations. 

Eq. ( [T^ ) represents a chain of coupled Langevin equations describing the evolution of 
interacting "particles". Even without solving the model, it is possible to realize that only 
one parameter suffices to describe the scaling behaviour of the Lyapunov exponent, the 
rescaled smallness parameter 

g = ^la\ (14) 

which can also be interpreted as the inverse effective diffusion constant. In fact, the factor 
7 in front of the deterministic forces can be eliminated by properly scaling the time units. 

In the next sub-section, we investigate numerically the validity of this prediction by 
computing the Lyapunov exponent for different values of a and 7 and thereby checking 
whether the data collapses onto a single curve. 

The evaluation of (5A(7) requires finding the invariant measure for the whole set of 
stochastic equations ([T^). This is still a difficult problem, since the deterministic forces do 
not follow from a potential and, therefore, the corresponding Fokker-Planck equation cannot 
be straightforwardly solved. In particular, it is interesting to notice that, although we are 
working in the limit of weakly coupled maps, the "particles" are not weakly interacting. 
This is the most serious difficulty towards a perturbative treatment of the problem. Nev- 
ertheless, we are going to see that this infinite set of stochastic equations can be effectively 
approximated by a single Langevin equation. 

Formally speaking, the probability distribution Pl{wi, ...,WL',t) satisfies the Fokker- 
Planck equation 

a p L Pj ( > 1 ^ pp. p 

^ - -g + *(u..«.u,._o/>4 + 2 S (15) 

where 

F(wi) = -27sinh(w,) (16) 
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is the single-particle force (see Eq. ([T^), = — 27(e'^*+i —e~'^''^) is the coupling 

term and T>i j are the diffusion coefficients as from ([T3|) . 

Let us now introduce the single-particle distribution P{wi) by integrating over all the 
other L — 1 degrees of freedom 



P{wi,t) = jY[dwj Pl(wi, ■■■,WL;t) . 



The corresponding equation can be directly derived from Eq. (|T 
BP 8 



Eq. (|T^) is not closed since it involves the unknown two-body distribution P2. The above is 
indeed the first of a hierarchy of equations involving probability distributions of increasing 
order. The simplest approximation to close the system consists in truncating the hierarchy 
at the lowest level by assuming a perfect factorization, P2{x,y) = P{x)P{y). This is the 
standard mean-field approach which leads to the closed Fokker-Planck equation 

BP B B'^P 

^._^{^HPK2.^|^, (18) 

since the symmetry of the distribution {P{w) = P{—w)) implies that the coupling terms 
cancel each other. 

The Fokker-Planck equation ( |T8|) corresponds to the single-particle Langevin equation 

w = F{w) + m, (19) 

where we have dropped the by now irrelevant spatial dependence. 

It is instructive to test the validity of the above approximation in the limit of large g 
values, when the forces can be linearized and an analytic solution can be obtained for the 
entire chain of Langeving equations. For the sake of readability, the discussion of this tech- 
nical problem is presented separately in the first Appendix. The main result of this analysis 
is that the approximation (p!9D is exact! The stationary probability distribution of Eq. (p!8D 
is equal to the projection of the many-particle distribution for the whole chain. However, 
one cannot expect that the same holds true also for finite g values, when nonlinearities come 
into play. 

Let us therefore discuss a possible line of thought to go beyond the approximation (|19|). It 
is first important to recognize the special nature of the noise ^, which keeps W(t) = J^i'i^iit) 
constantly equal to zero. Indeed, W is nothing but the logarithm of the ratio between the 
amplitude of the perturbation in the first and the last site. If periodic boundary conditions 
are assumed, as we do, this implies that W{t) = 0. A simple way to satisfy this constraint 
is to assume Wi = —Wi-i. By substituting this extreme assumption into Eq. (0), we obtain 
again a Langevin equation with the same force as in Eq. (|TBp but a different factor. It is 
therefore reasonable to conjecture that, in general, the effective force strength depends on 
g, thus writing (in rescaled time- units) 
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F{w) = —2a{g) sinh(w) 



(20) 



where a is a renormalization factor to be determined with some self-consistency argument. 
As we have not found any sensible way to do so, we limit ourselves to investigate numerically 
its possible dependence on g. 

The stationary distribution w is obtained by solving the Fokker-Planck equation 



P(w) = C exp[— a7/cr^ cosh(w)], (21) 
where the normalization constant C = 2Ko{a'~f /(t"^) is expressed in terms of the zeroth-order 



modified Bessel function 12,13 



/oo 
(it exp{— x cosh(t) ± vt}. 
-oo 



By substituting the definition of w {w = Ini?) in Eq. (|^), and assuming {wiWi^i) = (an 
hypothesis accurately confirmed by direct numerical simulations), we finally obtain 



As 7 is assumed to be small, we can safely retain only the first two terms of the series (P2D, 
obtaining 

which represents the (approximate) perturbative expression for the correction to the MLE 
in the limit of small coupling. 

First of all, it is instructive to investigate the limit (7 ^ 1, by using the asymptotic 
expression of the functions K,j{y). Since 

K,{y)^-\n{y/2) , K,{y) (24) 

y 

we find that 

W,) ~ ^ (25) 

This equation represents a relevant improvement over the previous results. First of all, it is 
in agreement with numerical simulations which do not give evidence of a In | ln(yf| correction 
in the numerator, as instead predicted by the statistical-mechanics treatment based on the 
analogy with directed polymers 0. 

A second and more important remark concerns the dependence on the "noise" strength 
which is explicitely determined. Previously it was only clear that the correction to the MLE 
must vanish if there is no multifractality (no multiplier fluctuations) but the dependence on 
cr was not known. 
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B. Numerical results 



The theoretical analysis performed in the previous section is mainly based on a perturba- 
tive approach. Moreover, it involves a nontrivial transformation of a set of coupled Langevin 
equations to a single effective Langevin equation in a limit where the coupling is not negligi- 
ble. Therefore, a comparison of the theoretical predictions with direct numerical simulations 
is worth especially to check the validity of the dynamical mean-field approximation that is 
behind this last step. 

We have decided to test the thereotical results by using two different probability densities: 

A) uniform distribution of multipliers rrii within the interval [e°(l — Ai/2),e'^(l + Ai/2)]; 

B) uniform distribution of the logarithms of the within [a — A2/2, a + A2/2]. The 
corresponding Lyapunov exponents in the uncoupled limit (e = 0) are 

-^<«)^"— ir{(-t)-(-t)-(-t)K-t)} 

As(0) = a, 

while the variances of Inmj(t) are 

9 1 / A2\ , 2 /! + Ai/2\ , , 



We start from testing the predictions for the shape of the probability distribution of w. Two 
meaningful examples are reported in Fig. |I|, where the open circles refer to the histograms, 
while the solid line is the theoretical result (Eq. (|2l| ) with a set equal to 1). Let us first 
comment about the qualitative shape of the distribution. In the limit of large g, the noise 
is almost negligible and therefore, the phase point is expected not to deviate significantly 
from the stable fixed point w = 0. It is therefore possible to linearize the equation, finding 
a Gaussian distribution. This is precisely the message contained in Fig. |l]a, which refers 
to g = 2.4. In the limit of small g instead, it is the force that can be neglected except 
when the deviations are large. Since the attracting force grows very rapidly (exponentially), 
it makes sense to replace the corresponding potential with a fiat well with infinitely high 
barriers placed where the deterministic force is of the same order as the stochastic one. In 
this picture, one expect that the probability distribution is just a uniform distribution in a 
finite interval (this is the kind of argument introduced in |^ to predict the scaling behaviour 
in this regime). This scenario can be qualitatively recognized in Fig. P^, which refers to 
g = 0.021. 

Next, let us comment about the quantitative agreement between the theoretical expec- 
tations and the numerical findings. In Fig. Qa, there is an almost perfect agreement. This 
is a first encouraging result, since it indicates that even for not too large a value of g, the 
effective force strength remains equal to 27 with no renormalization. Some deviations are 
instead observed in Fig. |l]b, where the theoretical curve is slightly more peaked, suggesting 
that the effective force is smaller than expected. A simple way to determine a is by fitting 
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the numerical data. In the first case we find that a = 0.97 confirming the first quahtative 
impression; in the second case a = 0.75 indicating stronger but not significative deviations 
from the hnear regime. The curve corresponding to the fitted value of a is reported in 
Fig. |l]b as a dashed line. The residual, small, deviations with respect to the numerical data 
indicate that the reduction of the original model to a single stochastic equation with the 
effective force (^Oj) is indeed a strategy that is worth being pursued. 

All the other cases that we have tested reveal the same scenario. Altogether, we can 
summarize stating that the major discrepancy between numerics and theoretical prediction 
is contained in the renormalization factor a which is assumed to be equal to 1 in our 
treatment but turns out to depend on g; nevertheless it is never smaller than 0.7. 

Since the aim of the present paper is to study the corrections to the MLE induced by 
the spatial coupling, let us now discuss this issue. In order to assess the quality of our 
theoretical predictions, we performed numerical iterations of Eq. (Q) computing the MLE 
with the well known algorithm . Simulations have been carried out on 500-site lattices. 



imposing periodic boundary conditions. In every simulation the first 500 iterations have 
been discarded to avoid any bias effect due to initial conditions. Tests made with different 
lattice lengths indicate that finite-size effects are always much smaller than the deviation 
from the theoretical predictions. 

The first nice result is provided by Fig. ^a. The data is plotted in order to emphasize the 
existence of only one relevant parameter, g. Indeed, the good data collapse (all data align 
along the same curve irrespective of the value of 7, cr or the type of probability distribution) 
represents a further confirmation of our theoretical analysis. Moreover, the nice agreement 
of the numerical data with the theoretical expression p3| ) (see the solid line in Fig. ^) over a 
wide range of values of the effective coupling testifies to the accuracy of the approximations 
introduced in the first part of this section. 

However, there is a better way to emphasize the differences between theory and numerics. 
In fact, the limit (7 ^ 00 corresponds to negligible noise, i.e. to a regime where the MLE 
Lyapunov exponent is unaffected by the presence of spatial coupling as shown in It 



is therefore convenient to look at the behaviour of the whole deviation of the Lyapunov 
exponent A A = 5A — 27, which again exhibits the same scaling behaviour as seen by 
dividing this expression by and using Eq. (^) 

AA_ //A(2j) \ 

The data plotted this way are reported in Fig. ^d. We clearly see that the trivial correction 
term —2'j cancels almost exactly the growth exhibited by SA for large (^-values allowing 
for a more stringent test of the theoretical prediction. We can now see that the absolute 
difference is not larger than 0.05 and it could be greatly reduced by suitably shifting the 
theoretical curve (i.e. by rescaling g by approximately a factor 2) as shown by looking at 
the dashed curve. However, in the absence of theoretical arguments, this observation cannot 
be considered more than a hint for future considerations. 



IV. THE GENERAL CASE 
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A. Theory 



In this section, we account for sign fluctuations as well. In order to keep the theoretical 
treatment as simple as possible, we shall assume that the sign is a (5-correlated stochastic 
process independent of the modulus, so that the average factorizes as 

{mi) = {p-q){\mi\), 

where p (q) is the probability that rrii is positive (negative). The major difference with 
respect to the previous case is that the ratio Ri can also assume negative as well as positive 
values. It is therefore more convenient to work with instead of introducing its logarithm 
which would require introducing absolute values and thus two different variables to account 
for the dynamics in the positive as well as in the negative range of Ri values. 

We have learned in the previous section that neglecting the coupling with the neighbour- 
ing sites provides a good approximation of the probability distribution and thereby of the 
MLE. Let us therefore neglect the dependence on Ri+i and in Eq. (|]). As a result, we 
obtain the one dimensional mapping 

where, for the sake of simplicity, we have removed the now irrelevant site dependence. In the 
continuous-time approximation, the above would be exactly the equation that has allowed an 
approximate analytical treatment in the case of positive signs. Mapping ( ^9]) possesses two 
remarkable symmetry properties: the evolution is invariant under the transformation R — »• 
1/R, since the stochastic process fi turns out to be invariant under the same transformation 
;U l//i (it is sufficient to look at its definition). This is the same symmetry as that one 
discovered in the previous section when we have found that the potential V{w) is an even 
function of w. As a consequence, we can restrict our analysis to the interval [—1, 1]. 

The second symmetry has much more serious implications. We can see that the mapping 
(^) is also invariant under time-reversal. More precisely, if we express R{t) as a function 
of R{t + 1) we find the same functional form of mapping (^), provided that the changes 
of variables S : R{t + 1) — > —R(t + l)//i(t) and fi 1/fi are introduced as well. As 
the transformation 5* is an involution (i.e., 5^ = Id), we can state that mapping ( ^^ is 
invariant under time-reversal. Therefore, we seem to be in the presence of a contradiction, 
as this property holds also for strictly positive multipliers, when it is clear that there is 
an attractor (the point w = 0, i.e. R = 1), since time-reversibility hints at a lack of 
attractors! Indeed, there is no contradiction, since time-reversal simmetry is broken in the 
case of strictly positive multipliers. In fact, invariance of the mapping implies only that a 
given solution can be neither stable nor unstable, if it is itself invariant under the involution 
S. However, this is not the case for the fixed point i? = 1 (in the absence of noise, i.e. for 
/i = 1) which is mapped by S onto R = —1, so that we can only conclude that if i? = 1 is 
stable, then R = —1 must be unstable (as it is indeed the case). In other words, positive 
values of R are characterized by a contracting dynamics towards R = 1, while negative 
values depart from —1. If the multipliers are strictly positive, negative i?- values cannot be 
asymptotically observed as they lie in the repelling part of the phase-space and the previous 
treatment in terms of a Langeving equation with an attracting force makes perfectly sense. 
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On the other hand, if the multiphers can assume both signs, the dynamical rule allows 
visiting interchangeably the positive as well as the negative region. In principle, it is still 
possible to have, on the average, a global contraction provided that a longer time is spent 
in the positive region. Actually, this is the assumption more or less implicitely made in 



Ref. |]T6[, where it was conjectured that no qualitative changes are expected when positive 
and negative multipliers come into play except for the degenerate case p = q = 1/2. We see 
below that even if the scaling behaviour in the limit of vanishing coupling is unaffected, this 
is not true and it indeed requires introducing a different scaling parameter. 

The most effective way we have found to analyse mapping (^) is by exploiting another 
property: the possibility to transfer the change of sign of /i to 7. With this trick, the change 
of sign in Eq. (|29| ) can be effectively treated perturbatively being 7 a small parameter. More 
precisely, if happens to be negative, we can assume it to remain nevertheless positive 
and perform the next iteration with —R{t + 1). It can then be seen that the resulting 
expression is the same as the original one after changing the sign of 7 and of + 1). Now, 
irrespective of the sign of —nit + 1), we assume it to be positive and transfer its sign to the 
next value of 7. In other words, we can iterate the mapping 

where the sign of 7(t) is that of 11^=1 /^(s)- We can immediately see that even if /i is on the 
average more positive than negative (or vice versa), the sign of 7 has no preference, since 
it simply depends on the parity of the number of sign changes. It is because of this reason 
that fluctuating multipliers are qualitatively different from strictly positive ones: even an 
asymmetry in the signs (a preference, say, for the positive values) implies that the unstable 
and stable region (positive and negative values of R in the initial representation) are equally 
visited. 

The dichotomic structure of the noise 7(t) allows expressing the stochastic map as the 
sum of a net drift plus a zero-average fluctuating term. Indeed, by calling -F+(-R) and F_[K) 
the l.h.s. of Eq. (^) whenever 7 is positive, respectively, negative, we can write 

R{t + 1) = \[F^m)) + F.{R{t))^ + ^{F+(i?(t)) - F_(i?(t))} (31) 
where 5{t) is again a dichotomic noise with entries equal to ±1. More specifically, we obtain 

In order to obtain an analytic expression for the probability density of i?, it is convenient 
to turn this equation into a continuous-time model. 

This is possible at the expense of assuming that the modulus fluctuations of ^ are small, 
i.e. be writing = 1 + z/, and by then expanding the r.h.s. in (0). By retaining terms up 
to second order, we obtain 

R = -f^R{R^ - 1) + (z/ + aV2)i? + (5(t)z/7(l - i?^), (33) 

where is the variance of u. 
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Before going on, it is important to notice that the above equation must be interpreted in 



the Ito sense as it arises from a discrete time stochastic process with 5-correlated noise |T7 
(this problem did not arise in the previous section, since we were deahng with an additive 
noise). Moreover, it is instructive to notice that the drift term in Eq. (|33|) is purely induced 
by noise: it arises from the inhomogeneity of the stochastic process. This represents a direct 
confirmation that contraction and expansions processes tend to compensate each other as 
already anticipated in the beginning of this Section. 

If the time variable is rescaled by a factor a^, it is immediately recognized that the 
dynamics of R depends on just one parameter 

G = =, (34) 
cr 

which is again a ratio between coupling strength and multiplier fluctuations. However, there 
is an important difference with the parameter (7 = 7/cr^ introduced in the previous case, as 
it is seen by noticing that in the small a limit, the equahty 



cr 



V2a (35) 



holds. Apart from the irrelevant numerical factor, it turns out that, in the general case, the 
r.m.s. rather than the standard deviation enters as a measure of multiplier fluctuations. 

The Fokker-Planck equation corresponding to the Langevin process (|33D in rescaled time 
units reads as 

dP d 1 92 

where 

A{R) = -G^R{1 - R^) + R/2 (37) 

is the drift term, while 

B{R) = R^ + G^{l-R^f (38) 
is the diffusion coefficient. Since 4:A{R) = dB/dR, the stationary solution is 

P,«) . ^-(g) (39) 

/G2(1 - i?2)2 + R2 



where N{G) is the normalization constant discussed in Appendix B. 

An expression for the Lyapunov exponent can be obtained from Eq. (|^), by integrating 
over the above determined probability distribution. Unfortunately, there is a crucial differ- 
ence with the previous case: we cannot simply expand the logarithm in powers of 7, since 
this leads to computing the first moment of P{R) which is already a diverging quantity 
{P{R) decays to zero as slowly as l/R"^). Obviously, this is only a numerical artifact: the 
average of the logarithm itself is still well defined and has a finite value. 

Nevertheless, this is an indication that we must be much more cautious in performing 
power expansions. In particular, this difficulty prevents obtaining a general analytical ex- 
pression analogous to that one obtained in the previous section in terms of modified Bessel 
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functions. In this case, even obtaining an expression in the hmiting case of small G requires 
a rather laborious work. In Appendix B we show that one can eventually show that the 
non-trivial deviation with respect to the uncoupled limit is given by 

Therefore, we see that also in the general case of positive /negative signs, the leading depen- 
dence on e is of the type 1/ Ine, as numerically observed. What is different is the dependence 
on the multiplier fluctuations as testified by the presence of the parameter G rather than g. 



B. Numerical results 



The first meaningful test of the analytical approach devised in the previous sub-section 
concerns the probability distribution P{R)- In Fig. ^we report the outcome of a numerical 
experiment in doubly- logairthmic scales (see the full dots). This allows seeing a crossover 
from an initial decay as 1/Rto the asymptotic decay which represents the first qual- 

itative confirmation of the thereotical predictions. However, the agreement with expression 
(p9|) (represented by the dashed line) is more than just qualitative. In fact, besides notic- 
ing the almost perfect overlap, one should also remember that the only parameter entering 
Eq. (^), i.e. G, has not been fitted, but independently determined from the fluctuations of 
the local multipliers. As a last remark, we would like to point out that the good agreement 
is not totally obvious a priori at least for the reason that the reduction from a set of coupled 
stochastic equations to a single equation is not completely under control. 

Moreover, it is instructive to compare the shape of this distribution with the results 
predicted by the theory for strictly positive multipliers. By expressing the probability density 
of Eq. (^Tj) in terms of R, we find an exponential tail {P{R) — exp{—gR)/R). The power 
law observed in Fig. ^ is also, therefore, an evidence of a clear difference between the two 
regimes. 

Finally, let us look at the deviations of the MLE plotted versus the scaling parameter 
G. The data reported in Fig. ^ have been obtained for different noise amplitudes and either 
e = 10~^ (diamonds) or e = 10~^ (circles). It is clearly seen that, SAlnG/a"^ is constant, 
independently of the value of G. This confirms the scaling behaviour predicted by Eq. (^OD. 
A more quantitative check can be made by comparing the actual value of 6 A In G/a^ (about 
1.1 ^ 1.2 in direct simulations) with the theoretical prediction (1.5) . We believe that 
the deviation is to be ascribed to the approximation made in reducing the set of coupled 
stochastic equations to a single Langevin equation. 



V. CONCLUSIONS AND PERSPECTIVES 



In this paper we have developed a theoretical method that is able to explain not only 
the scaling behaviour of the maximum Lyapunov exponent for a CML in the small coupling 
limit, but provides a quantitative estimation of its deviations from the uncoupled case. This 
was still lacking even in the relatively simple case of strictly positive multipliers. However, 
we have gone further, developing a treatment also for multipliers with random signs. In 
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both cases we have found that the correction to the maximum Lyapunov exponent induced 
by local interactions actually depends on a single scaling parameter which is nothing but the 
coupling strength rescaled by the "amplitude" of multiplier fluctuations. However, the scal- 
ing parameter is significantly different in the two cases: for strictly positive multipliers, the 
fluctuation "amplitude" is the mean square deviation (see the definition of g - Eq. (p!^), 
while in the case of random signs, the "amplitude" is the r.m.s. deviation (see the definition 
ofG-Eq. (0)). 

Important differences can also be detected in the probability distribution of the local 
ratios Ri of the perturbation amplitude in two adjacent sites: in the case of fluctuating signs 
there are long tails characterized by a power law decay. 

Among the still open problems, there is certainly the exigence of a more rigorous proce- 
dure to solve the set of coupled Langeving equations. In fact, while the derivation of the set 
of coupled stochastic equations is the result of a well controlled perturbative approach, its 
reduction to a single equation is based on a mean-field approximation whose validity cannot 
be controlled a priori but only checked a posteriori. 

Finally, we want to mention the possibility of extending this approach to the case of 
weakly coupled attractors, where time is continuous from the very beginning. This is cer- 
tainly the most stimulating perspective that is also supported by the preliminary observation 
that the Langevin equation (O) is obtaind also in the case of two weakly coupled differential 



equations [11 
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APPENDIX A: LINEAR LIMIT 

This Appendix is devoted to the analysis of Eq. (|T2|) in the linear limit, 

Wi = -'j{wi+i - 2wi + Wi^i) + ipi{t) - ^i-iit), (Al) 

where we have introduced ipk = ln{mfc(t)}. This is apparently a discretized Edwards- 
Wilkinson equation [|1^], but the spatial structure of the noise prevents the onset of any 
roughening phenomenon (as commented in the main body of the paper). 

To solve this equation, it is convenient to perform a spatial Fourier transform, since it 
leads to a set of uncoupled equations, 

w{k, t) = -27(1 - cos{k))w{k, t) + (1 - e''')^{k, t), (A2) 

where w{k,t) is a complex number that can be decomposed into a real and imaginary part 
{w{k,t) = x{k,t) + i y{k,t)), which satisfy the same equation 

x{k, t) = -27(1 - cos(A;))a;(A;, t) + r]{k, t), (A3) 

where the noise term 77 is 5-correlated, 
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{ri{t)ri{t')) = 2a^ (1 - cos{k))5{t - t'). 



(A4) 



Accordingly, all Fourier modes obey the same Gaussian distribution function, 

P{.}~..p{-'-^^). (A5) 

The probability distribution of Wi on a single site is easily obtained by summing the inde- 
pendent distributions corresponding to all modes. As a result, the distribution of Wi is also 
Gaussian and its variance is cr^, as if we had neglected the spatial coupling in Eq. ([Al|). 



APPENDIX B: LYAPUNOV CORRECTION 



In this Appendix we determine the nontrivial part of the leading correction to the MLE 
in the general case. We start by computing the normalization constant. It is convenient to 
exploit the invariance of P{R) under parity change and the transformation i? — > to 
express the normalization condition as 

1 = 4/ dRP(R) = 4N(G) / . (Bl) 

•^0 ^0 ^R2 + ^2(1 - i?2)2 

Since an explicit analytical expression for the above integral does not exist, we shall limit 
ourselves to studying the small-G limit. One cannot simply expand the denominator, as it 
gives rise to a non-integrable singularity in = 0. It is, instead, convenient to introduce 
the variable x = R/G. Afterwards, one can expand the integrand in powers of G without 
encountering undesired divergences. By retaining the leading terms, we find 

1-4 f^''^iL=^41nfl/G'') (B2) 



N Jo 



X 



From Eqs. (||,|), it turns out that the estimation of 6 A requires computing the mean value 
of L(i?i,i?2) = ln|l + 7/?i + 7/i?2|, i.e. 

/OO roo 
/ dRidR2P{Ri)P{R2)L{Ri,R2) (B3) 

Thanks to the equality L(i?i,i?2) = L{1/ R2,l/ Ri) and to the invariance of P{R) under 
the transformation i? — > 1/i?, we can write the Lyapunov correction as the sum of three 
different contributions, namely 

6A = 61 + 62 + 63 (B4) 
= J'^ dRidR2P{Ri)P{R2)[L{R^, I/R2) + 2L(i?i, R2) + Lil/R,, i?2)| 

with an obvious meaning of the new symbols. 

In analogy with the computation of the normalization constant, we introduce the vari- 
ables X = Ri/G and y = R2/G. As a consequence, the expressions for the three contributions 
write as 
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* = r r dx<;,^^te£i±2£^_ (B5) 

J-l/G J-l/G + (1 - G2a;2)2y^2 + (1 _ Q2y2y 

= r r , ■°l^ + ^'^^t^^<'^^'l (B6) 

ii/G ./a;2 + (1 - G^x^yJy^ + (1 - G2y2)2 



= , ■n|l + 0V(g^)+7/(Gy)l ,B7) 

7-1/G J-l/G ,/x2 + (1 - G^x^yJy^ + (1 - GV)^ 



The inequalities 7 <^ G ^ 1 imply that the contributions proportional to 7G in the argu- 
ments of the logarithms can be neglected so that 61 is negligible altogether. Moreover, G 
can be always neglected in the denominators, so that the leading contribution to the MLE 
can be determined by just estimating the two integrals 

/•i/G ln|l + CT/y| 
6. = N' r rdW^^'^^ (B9) 

J-l/G J-l/G ^ ^/TT^yrTF 

where we have re-introduced the parameter a for later convenience. We start discussing ^2; 
it cannot be computed by expanding the logarithm as this leads to an unphysical divergence. 
It is, instead, helpful to split this contribution into two parts 

62 = 62 + 62 (BIO) 

/■i/G Inlv + ffl /"i/G , Inlyl 

= 2N dy , ,' - 2N / dy- 



-i/G VTTy^ J-l/G VI + y'^ 

The first integral can be estimated by introducing the variable w = y + 'a and thereby 
expanding the denominator in powers of a. By retaining terms up to the second order, we 
find that 62 can be written as 

6'2=2N dw\^=J=J= + a——Ajrin:, + —--, !j-^{2w^-l)\ (Bll) 



(1 + ^2)3/2 2 (1+^2)5/2' 



By expanding the zeroth-order term around the integral boundaries in powers of a, we find 
that it is equal to —62 plus corrections of the order 7^ InG. A contribution of the same order 
is obtained also by integrating the linear term in a. However, the leading contribution to 
the MLE comes from the second-order term which, in the small-G limit can be written as 

S2 = 2Na^l ^-(Y^^(2-^-l)- (B12) 



The integral can be analytically solved and turns out to be equal to 1, so that [|T9 

The determination of ^3 in principle requires even more cumbersome calculations, as it 
involves a double integral. However, formally deriving the expression for 63 /N"^ with respect 
to G, we find that, up negligible corrections. 
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GJ-i/G VT+f 

As the integral in this expression is exactly the same involved in the definition of 62 (see 
Eq. (pHD), we can write 

Upon substituting the expression for 62 (see Eq. ( |B13|) , the above equations can be rewritten 
as 

which, after integration, yields 
In conclusion, we find that 

where we have preferred to introduce the explicit dependence on the physical parameter a 
rather than a. 
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FIGURES 



FIG. 1. Probability distribution of w for two different values of g: 2.4 (a) and 0.021 (b). In both 
cases, circles refer to the numerical histograms, obtained by iterating Eq. (||) for 5 x 10'' time-steps 
(discarding the first 10'^ iterations) on a lattice of L = 300 sites. The solid curves correspond to 
the analytic formula (Eq. (|2T| ) with a = 1). The dashed curve in (b) is obtained by fitting a which 
is estimated to be equal to 0.75. 

FIG. 2. The Lyapunov exponent versus the scaling parameter g in the case of strictly positive 
multipliers. The data is obtained by varying e, a and for both choices of the probability distribution 
of the local multipliers (the cases A and B discussed in the text). In a), the shift 6X defined in 
Eq. (|^ is reported, while in b) the total shift A A (see Eq. (p^)) is plotted. The solid curves 
correspond to the analytic expression. The dashed line in b) is the analytical curve arbitrarily 
shifted to show that a "renormalization" of the scaling parameter could account for the remaining 
discrepancy with numerical data. 

FIG. 3. Log-log plot of the probability distribution P{R) to highlight the power-law behaviour. 
Circles, triangles and diamonds refer to G = ■ 10~^,\/6 ■ W~^,VQ- 10~^, respectively. The 
simulation details are as in Fig. |l[ The various curves represent the analytical results as from 
Eq. (H). 

FIG. 4. The Lyapunov exponent versus the scaling parameter G in the case of fluctuating 
multipliers. The Lyapunov correction 5 A is normalized so as to emphasize the 1/| In G\ dependence. 
Circles correspond to e = 10~^ while diamonds to e = 10~^. The straight solid line represents the 
theoretical result (|^). 
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